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Abstract. 

It is shown that the magnitude of tensor polarization of the deuteron beam, which 
arises owing to the spin dichroism effect, depends appreciably on the angular width 
of the detector that registers the deuterons transmitted through the target. Even 
when the angular width of the detector is much smaller than the mean square angle 
of multiple Coulomb scattering, the beam's tensor polarization depends noticeably on 
rescattering. When the angular width of the detector is much larger than the mean 
square angle of multiple Coulomb scattering (as well as than the characteristic angle of 
elastic nuclear scattering), tensor polarization is determined only by the total reaction 
cross sections for deuteron-nucleus interaction, and elastic scattering processes make 
no contribution to tensor polarization. 
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1. Introduction 

The quasi-optical macroscopic quantum phenomenon of birefringence, first described 
in [U [2] , combines the effects of particle spin rotation (oscillation) and spin dichroism 
(i.e., the dependence of the absorption coefficient on the particle spin state) and occurs 
when high-energy particles with spin S > 1 traverse matter. It is analogous to the 
phenomenon known in optics as birefringence of light in optically anisotropic media due 
to the dependence of the refractive index on the polarization state of light, i.e., on the 
photon spin state. In contrast to light, whose wavelength largely exceeds the interatomic 
distance in matter, the de Broglie wavelength for a fast particle is much shorter than 
this distance. According to [U [2j [3], however, in this case one can also introduce for 
particles a spin-dependent index of refraction. Under such conditions, the birefringence 
effect for deuterons takes place even in homogeneous, isotropic matter and results from 
the intrinsic anisotropy of particles with spin S > 1. 

An important implication of the spin dichroism effect is tensor polarization arising 
in the initially unpolarized beam after it has passed through an unpolarized target [TJ [2] . 

In 2005, the spin dichroism effect was first experimentally observed for 5-20 MeV 
deuterons passing through a carbon target jUGEE]. In 2007, tensor polarization arising 
in a nuclotron extracted beam of initially unpolarized deuterons with a momentum of 
5GeV/c was measured [HE]- 

The magnitudes of spin dichroism and tensor polarization of the transmitted 
deuterons are proportional to the difference between the total interaction cross sections 
of the particle with the nucleus for the states with magnetic quantum numbers M = 
and M = ±1 PEE]. 

According to p!J|2l|3jlIJ|5j|6j[7JIE], this effect can be used not only for studying the 
total spin-dependent interaction cross section, but also for obtaining tensor-polarized 
deuteron beams. Besides, this effect must be taken into account when interpreting the 
results of planned experiments [9] for measuring the electric dipole moment of deuterons 

mm- 

In theoretical description of this phenomenon, one should bear in mind that a 
deuteron colliding with a nucleus undergoes two interactions: the Coulomb and nuclear 
ones. The Coulomb interaction not only causes multiple Coulomb scattering, but in some 
cases, it can also affect significantly the magnitude of tensor polarization of the deuterons 
that have passed through the target. It has been shown in [UJ, for example, that 
the Coulomb-nuclear interference gives a clue to the explanation of the experimentally 
observed fact that tensor polarization reverses sign as the energy of deuterons changes 
from 5 to 20 MeV. 

Multiple scattering of deuterons in matter and the scattering-angle dependence of 
the Coulomb-nuclear interference account for the necessity to use the kinetic equation 
for the spin density matrix in order to describe the spin-state dynamics of the deuteron 
passing through matter [12] . 

In the present paper, it is shown that tensor polarization of the deuterons escaping 
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from the target in the direction of the in-coming primary beam depends on the 
collimator's angular width of the detector that registers the transmitted deuteron 
beam. It is demonstrated that when the angular width of the detector's collimator 
is much larger than the characteristic angles of deuteron scattering in the target, tensor 
polarization of deuterons is determined only by the total cross sections of inelastic 
processes in the interaction between the deuterons and the nuclei. 

This paper is arranged as follows. Section [2] gives the kinetic equation of the 
spin density matrix and the spin-dependent nuclear scattering amplitude. Section 13.11 
considers deuteron spin dichroism in a thin target, where only single scattering is 
essential. Section 13.21 discusses the effect that multiple scattering makes on the 
magnitude of tensor polarization of the deuteron beam moving in matter. A brief 
summary of the obtained results is given in Conclusion. 



2. Kinetic equation for the spin density matrix 

In order to describe the behavior of a beam of spin particles in matter, let us introduce 
the spin density matrix p of the system "incident particle + target". The spin density 
matrix pd of particles (deuterons, protons) moving in matter is defined by the expression 
Pd = Trx p, where Trx means summation of the diagonal elements of the matrix over all 
states of the target. The equation for time evolution of the spin density matrix pd has 
the form [13] : 

dp d ir- A i , fdpA 



dt h 



+ 



dt 



set 

where the Hamiltonian H describes the interaction between the particles and 
macroscopic external fields, and the collision term (dpd/dt) sct describes the density 
matrix evolution due to scattering by the target nuclei (atoms). 

The explicit form of (<9pd/<9t) sct can be obtained when considering scattering of the 
incident particle by the bound target particle, provided that the following inequality 
holds: Wo/ek <C 1. Here W is the binding energy of the scatterer in the target, and e k 
is the kinetic energy of the incident particle. In this case, impulse approximation applies 
to calculating the scattering process [H]. We shall further consider the case when the 
external electric and magnetic fields are absent. This enables one to write the below 
integro-differential representation of the kinetic equation for the density matrix that 
describes the behavior of a polarized particle beam as it passes through matter whose 
nuclei (atoms) in the general case are polarized [12] : 

^dT^ = W r NTlT {Hl^,t)-/$,t)P+{k,ty 

+^Tr T J d 3 k' 6 (^e k - e k , - ^ F(k,k')p(P } t)F + (k' } k) } (2.1) 

where F(k, k') is the amplitude of particle scattering by the target nuclei, N is the 
number of target particles in the unit volume, M is the mass of the scatterer, M r is 
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= N Tr T 

dt m 

k 



the reduced mass of the incident particle and the nucleus, is the kinetic energy of 
the incident particle, k and k' are the wave vectors of the beam particle before and 
after scattering, respectively. The density matrix p{k) is the direct product of the spin 
density matrix p d {k) of the incident particle and the spin density matrix p^{k) of the 

target nuclei: p(k) = p d (k; S d ) <8> Pt(St), where S d and St are the spin operators of the 
beam particle and the scatterer, respectively. 

Kinetic equation (12.11) applies to such values of the transmitted momentum \q\ = 
\k' — k\ that are governed by the inequality q ^> r _1 , where r is the correlation radius 
of the medium [12J. 

Equation (12. ip simplifies when a particle (proton, deuteron, antiproton) passes 
through a target whose nuclei have a mass much greater than the mass of the incoming 
particle. In this case, we can neglect the effect of the incident particle's energy loss 
through scattering. So we can neglect the recoil energy q 2 /2M in the 5-function. As a 
result, one obtains a simple kinetic equation describing the time and spin evolution of 
the incident particle as it passes through the target [121 H3] : 

dp d (k,t) 27ri Ar ^ [P(k,k)p(k,t)-p(k,t)F + (k,k) 

+N- Tr T / d%F(£, k')p(k' } t)F + (k' } k), (2.2) 
m J 

where \k\ = \k'\ and m is the mass of the incident particle. 

The first term on the right-hand side of (12.21) can be represented as follows: 

F(0)p(M)-p(M)F + (0) 

= ~ (^(0) + F + (0)) ,p(M)] + {- (F(0) -F+(0)) J(k,t)}, (2.3) 

where [, ] is the commutator, { , } is the anticommutator. 

The part proportional to the commutator leads to rotation of the polarization vector 
due to elastic coherent scattering [13] (as a result of the refraction effect [U 121 E]); the 
anticommutator describes the reduction in the intensity of the transmitted beam. 

The last term in (12. 2p determines the effect of incoherent scattering on the change 
of pd (in the general case, the effect of single and multiple scattering). 

For the sake of concreteness, let us consider the process of deuteron passage 
through the target with spinless nuclei. In this case, the density matrix p(k), as 
well as the amplitude F(k,k'), contains only spin variables of the scattered beam. 
For the amplitude F(k,k'), we shall introduce the notation f(k,k') = F(k,k r ), where 

As a result, (refl.2) can be written as follows: 



(/(*, k) + ;+ 0, k)),p d (k, z)\ + jN{(f(k, k) - r(k, k)),p d (k, *)} 



dpd{k,z) = m N 
dz k 

+N J dQpf(k t k')p d (k', z)f + (k', k), (2.4) 

where z = vt (v is the particle velocity) is the distance travelled by the incident particle 
in matter. Hereinafter, the subscript d of the density matrix will be dropped. 
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In the case of deuterons (particles with spin 1), the polarization state is 

characterized by the polarization vector P(k) = Tt p(k)S and the quadrupolarization 
tensor (tensor of rank 2) Q, whose components are defined as Qi k (k) = Ti p(k)Q ik , 
where the operators Qi k can be represented in the form: Q ik = § (^SiS k + S k Si — ^8 ik Tj. 
The spin density matrix p(k) can be written in the following general form: 

p(k) = l I(k)\ + \p(k)§+ l -Q %k (k)Q lkl (2.5) 
o z y 

where I(k) = Tr p(k), I is the 3x3 identity matrix in the spin space. Note that alongside 
with the quantities P and Q, normalized spin characteristics of the beam are also used 
to describe the beam polarization [T5J HE] . 

For particles with spin 1 that are scattered by an unpolarized nucleus, the 
amplitude f(k,k') can be expressed in terms of the deuteron spin operator S, the 
quadrupolarization tensor Qi k , and the combination of vectors k and k'\ 

f(k, k') =AJ + B(Sv) + CxQikfiifik + C 2 Q ik p,xi^\k, (2.6) 

where A, B, Ci, and C 2 are the parameters depending on the scattering angle 9; 
v= [jfex k']/\[k x k'}\; (2= (k-k')/\k-k'\, and ft = (jfe + k')/\k + k'\. 

In view of (12.61) . for the zero-angle scattering amplitude f(k, k), we have 

f(k,k) = f (0) + W)(Sn) 2 } (2.7) 

where n = k/k is the unit vector in the direction of k, and the following notation is 
introduced: /„ = A - 2C X - 2C 2 , fi = 3C 2 - 

In the general case, fo and fi are complex functions, and according to the optical 
theorem, the imaginary parts of /o and f\ can be expressed in terms of the corresponding 
total cross sections, namely, 

lm/ (0) = A a t ° ot , ImA(O) = A [a^ - a° tot ] , (2.8) 

where cr° ot , cr to J are the total scattering cross sections for the initial spin state of the 
deuteron with magnetic quantum numbers M = and M = ±1, respectively (the 
quantization axis z is directed along n). 

Deuterons interact with the target nuclei via the Coulomb and nuclear interactions. 
The amplitude f(k,k') of deuteron scattering by the target nuclei in this case can be 
represented in the form [12j H7] : 

f(k, k') = /coul(fc, k') + / nu d,coul(fc, fc'), (2.9) 

where / cou i is the amplitude of Coulomb scattering of the deuteron by the nucleus in 
the absence of nuclear interaction, / nU ci,coui is the amplitude of scattering of Coulomb- 
distorted waves by the nuclear potential. 



3. The phenomenon of spin dichroism 



The phenomenon of deuteron spin dichroism (i.e., the dependence of the absorption 
coefficient on the spin state of the particle moving in matter) leads to the appearance 
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of tensor polarization in the initially unpolarized beam transmitted through the 
unpolarized target. 

The unpolarized deuteron beam can be considered as an incoherent mixture of three 
polarized beams having the same intensity and spin states with magnetic quantum 
numbers M = —1, 0, and +1, respectively (the quantization axis is directed parallel 
to the deuteron momentum). The total scattering cross section cr tot of deuterons by 
a nucleus depends on the polarization state of the incident particle; namely, in the 
interaction with an unpolarized nucleus, cx^t = 0"tot a tot- As a result, the partial 
intensities J (±1) and J (0) of the beam transmitted through the target appear to be 
unequal, and so the spin dichroism phenomenon can be characterized by the parameter 
n±i) _ r(o) 

D = - —, (3.1) 

J(±l) _|_ J(o) ' v 1 

where J (0) , J (±1) are the intensities of the deuteron beam transmitted through the 
unpolarized target for the case when the initial state of the deuteron is defined by 
quantum numbers M = and M = ±1, respectively. 

Spin dichroism leads to the appearance of tensor polarization in the transmitted 
beam. According to the definition (see e.g., [151 [E])> the diagonal components of the 
normalized tensor polarization can be written in the form: 

I /(+!) _|_ /(-!) _ 2/ (0) J (+1) + 7 (_1 ) — 2J ((,) 

Pxx — Pyy — — ^ j (0) 5 Pzz ~ j (+ i) j ( _i) j (0) ' 

where p xx = Q xx /I, p yy = Q yy /I, and p zz = Q zz /L 

In the case when the beam passes through the unpolarized target that is considered 
here, the intensities J (+1) = J (_1) . As a result, we have the following expression for p zz : 

J(±l) _ J(0) 

Pzz = 2 j , (3.3) 

where I is the intensity of the initially unpolarized beam, / = + + J (0) = 
2/(±D + jo). 

3.1. Tensor polarization of deuterons passing through a thin target 

Let us consider particle transmission through the target with thickness z, which is 
smaller than the deuteron's free path in matter, i.e., z < 1/Na, a is the total scattering 
cross section of the deuteron by the nucleus. In the first-order perturbation theory, the 
solution of kinetic equation ( 12. 4ft can be represented in the form: 

p{k, z) = p(k, 0) + jN [(/(£, k) + f + (k, k)),p(k, 0) 

+jN {(f(k, k) - f + (k, k)),p(k, 0)}z + Nf(k, k Q )p(0)f + (k o , k)z, (3.4) 

where p(k, 0) is the density matrix of the beam at entering the target, i.e., when z = 0. 
It describes the momentum distribution of the incoming particles with respect to the 
direction k . 
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In deriving ( 13. 4p . it was also assumed that the beam's initial angular distribution 
is much smaller than the characteristic angular width of the differential scattering cross 
section. * *. case, *. a.pHtude f(k, k>) h *. to Jdn- t ,f(k,k')p(k',0)f+(k',k) 

can be removed from the integral sign at point k' = k . As a result, 
we have J dQpf(k, k')p(k', 0)f + (k', k) ~ f(k,k )p(0)f + (k ,k), where p(0) = 

J dD,£,p(k' ,0) is the spin part of the beam's density matrix p(k,0). The term 

N J d%/(jfe, 0)/ + (#, fc) or JV/(jfe, k o )p(0)f + (k o , k) describes the contribution to 

the evolution of the density matrix due to single scattering of particles in the direction 
of vector k. 

Using the solution ( 13. 4ft and the explicit form of the spin structure of the amplitude 
f(k, k') [see (|2.6|) ]. one can find the dependence of the beam's intensity and polarization 
characteristics on the direction of particle scattering and on the distance z passed 
by the deuteron in matter. In a real experiment, the collimator of the detector has 
a finite angular width, and so the scattered particles are registered within a certain 
range of finite momenta. For this reason, the characteristics of the beam transmitted 
through the target should be studied in the range of solid angles AQ with respect to 
the initial direction of beam propagation (see figure [TJ). If the collimator has axial 
symmetry, then AQ is defined by the angular width of the detector's collimator 
(for further calculations, it is also assumed that i^det is much larger than the initial 
angular distribution of the beam). 




\ 



target 



Figure 1. Scheme of detecting the transmitted beam. 



We shall further consider the case when the unpolarized beam is incident onto the 
target. 

The number of particles X (0) , X (±1) with quantum numbers M = and 
M = ±1, which are registered by the detector with angular width AQ 



is determined by the diagonal matrix elements 



An 



d^p(k,z)\^ M=0 ) and 



An 



dVL^,p{k, z)\ty M= ±i) of the spin density matrix (the quantization axis is 



chosen along the direction n = k/k): 



l m (z) 



r(0) 
J 



1 - Na° tot z + N J dn (* u=0 \ (^\ \* M=0 )z 
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r(±l)/ 



r(±i) 

-'o 



1 - Nogz + N 



(3.5) 



where X = / dOJ(^, z); Iq 0) and J^ 1 ' are the partial values of the flow of the deuteron 

J An 

beam at z = for the spin states \^f M =o) and \^m=±i), respectively. For the case of 
unpolarized beams, considered below, Iq 0) = |/ and Iq 15 = |/o, where Jo is the deuteron 
flow at the target entrance; I = /dM(M) . The differential scattering cross section 
of the deuteron is (d<r/df2) sc . 

To find tensor polarization of the registered beam, let us make use of (13. 3p and 
(13. 5p . As a result, one can obtain the following expression for p zz : 



Pzz ^ ~^N(a 



2 - %~ tot 



that is, 



4n 



3 J ah 



/T , ^ 



dQ { (* M=±1 



/da N 

o) 



|**=±l) 



(3.6) 



Pz 



{o? - o°) z 



2 

— N 
3 



{of-O^)- / dfi (<* M=±1 | 
J An 



9 I*— 



(3.7) 



where <r r is the inelastic part of the total cross section, which includes all possible nuclear 
reactions, a e \ is the total cross section of elastic scattering. 

As is seen, taking account of the scattered particles leads to the presence in 
the tensor polarization of the term depending on the angular width of the detector's 

collimator. Particularly, if $ det — > 0, then / dVt \ — — ) — > 0. As a consequence [see 

J An \dil J sc 

(13. 6p ], tensor polarization depends on the difference between the total cross section for 
the deuterons in the initial spin state with M = and that with M = ±1: 

Pzzi^ct -+ 0) ~ ~N (<r& - a° tot ) z. (3.8) 

If i?det — ► TTj then all the scattered particles are registered. In this case, the 
second term is proportional to the difference between the total cross sections of elastic 
scattering: 



±i 

el 



(3.9) 
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As a result, according to (13. 7p . tensor polarization p zz is described by the following 
expression: 

jU^det = tt) ~ ~N (af - a?) 2. (3.10) 

The differential cross section of elastic scattering for fast particles scattered by the 
nucleus achieves the largest values in the range of angles 9 < 1/kR (R is the radius of 
action of forces) and decreases rapidly with growing 9 (9 3> 1/kR). For this reason, p zz 
achieves the value (I3.10P already when $ dct ^> 1/ kR. 

Thus, in the range of large values of i?<i et , tensor polarization only depends on the 
difference between the inelastic scattering cross sections for deuteron states with M = 
and M — ±1. It is quite understandable because in this case, there is no elastic- 
scattering-related losses in the flow of particles registered by the detector: the particles 
that have passed through the target without being scattered and those that have been 
elastically scattered equally get into the detector. 



3.2. Influence of multiple Coulomb scattering on tensor polarization of deuterons 
passing through matter 

The solution ( 13.41) applies to such target thicknesses z, for which multiple scattering in 
matter can be neglected, i.e., for z < 1/No~ cou \. 

To solve (12 .4p for the case z > l/Na cou \, let us take into account that in the 
considered energy range, <7 nuc i <C cr cou i. As a result, for the zeroth approximation we use 
the solution of kinetic equation (12.41) . where the collision term is only determined by the 
the Coulomb interaction between the incident particle and the nuclei of matter. The 
contribution coming to the evolution of the spin density matrix from nuclear scattering 
and the Coulomb-nuclear interference can be considered as a correction, which is true 
for the cases when z < l/No~ nVLC \. 



3.2.1. Solution of the Kinetic Equation. Using the optical theorem (12.81) . equation 
(12. 4p can be written as follows: 

a r o ~ / r* \ N 



Na° tot p(k,z) - - (o& - a° tot ) {(Sn) 2 , p(k, z)} 

^JVReA(O) [(Snf,p{k,z)] +N J dn p f{k,k')p(k',z)f + (k',k) (3.11) 



dz totFV ' ' 2 

2th 



where n is the unit vector in the direction of the momentum k. 

It follows from the form of the amplitude (I2.9P that the elastic part of the total 
cross section a'^ 1 can be presented as a sum of the terms describing the Coulomb 
scattering cross section, the nuclear cross section, and the contributions to the cross 
section coming from the interference between the Coulomb and nuclear interactions: 
(T^ 1 = cr° *J + o"° u *J coul + a° u *}. Let us introduce the following notations: cx^c 1 = 
a nucicoui + °iiug1 + ^"r'* 1 - Thus, the cross section a^ 1 can be presented in the form: 

^tot = ^coul + ^NC • I 3 - 12 ) 
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In view of ( I3.12p and the representation (12.91) for the scattering amplitude, kinetic 
equation (13.111) can be written in the form 



dp(fc,z) 
dz 



N 



-Na° coul p(k, z)-- (C ul - cr c ° oul ) {(Sn) 2 ,p(k, z)} 



2tt\ r -> -> 

+—NRefr\0) [(Sn) 2 ,p(k,z) 



+N / dQ p f coul (k,k')p(k',z)f+ ml (k',k) 



N 

-Na° NC p(k,z)~-(a 



±i 

NC 



a 



NC, 



{(Sn) 2 ,p(k,z)} 



+— iVRe/r cl (0) \(Sn) 2 ,p(k,z) 



+N J dQ p (f coul (k,k')p(k',z)f] 
+ /nuci,coui(&> k )p(k , z)fc OU \(k , k) 
+N ' 



nucl,coul 



(k', k) 



(3.13) 



dfig//nucl,coul(^, k )p(k , ^)/n UC l iCO ul(^ ; k). 

As stated above, the solution of (13 . 1 3[) can be presented as follows: 

p(k,z)=p^(k,z)+p^(k,z) + ... : (3.14) 

where p^°\k,z) is the zeroth approximation, which is the solution of kinetic equation 
( I3.13P in the case when only the Coulomb interaction between the deuteron beam and 
the nuclei of matter is taken into account; p^'(k, z) is the first-order perturbation theory 
correction, which includes the contribution of nuclear scattering to the evolution of the 
beam's polarization characteristics. 

The equation for p^(k, z) has the form: 



dp(°>(fc, 
dz 



N 



{(J. 



±i 



coul 

a(0) 



coul. 



-Na° oul p(°\k,z) 

+ ^ATRe/ 1 coul (0) {Sn) 2 , p m (k, z) 

k 1 



{(^) 2 ,p(°) 



(k,z)} 



(3.15) 

where the spin structure of f CO u\(k, k') in the general case has the form (12.61) . 

In the small-angle approximation, the terms in (12.61) that are proportional to B, 
Ci, and C 2 in the Coulomb amplitude, for k' ^ k lead to depolarization of the registered 
beam. The magnitude of this depolarization is determined by the quantity b 2 ; 8 2 z 
2 7 2 



9 



1 7 - 1 — 

- H , g is the gyro magnetic ratio, and 9 2 is the mean square angle 



2 7 7 

of Coulomb scattering per unit length) [3j [18]. For the considered path lengths, which 
are of the order of the nuclear collision length or smaller, the degree of depolarization for 
deuterons having energies in the range of hundreds of megaelectron-volts and moving 
in carbon is ~ 10 -3 . As is seen, deuteron depolarization is insignificant and will further 
be neglected. As a result, in the case of Coulomb interaction, it is sufficient to consider 
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the spin-independent part of scattering amplitude. Then (13.1 5p reads: 

' ? ' ^ = -Na° covl p^(k,z) + N f dQp a(k,k') 2 p (0) {k',z), (3.16) 



dz 

where a(k, k') denotes the spinless part of the amplitude f CO ui(k, k'). 

Expression (13.161) is an integro-differential equation. In the limit of small scattering 
angles, it can be solved by expanding the function p(°'(k') in terms of a small parameter 
q = k' — k<^ik (transmitted momentum) |19j : 

p^(k', z )^p^(k, z ) + ^ qx + ^ qy 

+ 2^T^' + dk~dk- y q * q y + 2~Wf q y + (3 - 17) 

where q x = gcos<p, q y = qsimp, and q^kQ with 9 being the scattering angle (the angle 
between vectors k and k'). 

Substitution of the expansion (13 . 1 Tf) into (13.161) and integration over the azimuthal 
angle <p gives 

= ?V°)(M) + ^AApM(k,z) + (3.18) 

8 2 d 2 

the operator A affects the transversal components of vector k: A = — — - + — — -, where 

an z x any 

k/k, "F = N / 9 2 ^±dn, and "F = N [ 6 4 ^±dtt, etc. 

J dil J dil 

If we confine ourselves to the first term on the right-hand side of (13 . 18[) , then integro- 
differential equation (I3.17P can be reduced to a differential equation. The limiting angle 
in such approximation is obtained from the condition 

16A,(°>(M) 

As a result, we have the equation for the spin density matrix p^°\k, z), which coincides 
in form with the equation for the spin density matrix describing multiple scattering of 
spinless particles: 



n 



^ = ^V>(M). (3,0) 

For the initial condition p(°'(k,z = 0) = po5(n x )5(n y ) and infinite media, the solution 
of this equation can be presented in the form: 

p(°\k,z) = p g(k,z), (3.21) 

where po is the spin part of the density matrix [see (12. 5p 

PO = — -^ol + 2^°*^ gQoikQik, (3.22) 

and the function 
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The unit vector n = k/k is counted off with respect to the direction of n (further, we 
shall align the z-axis with vector n ). In the small-angle approximation, (n — n ) 2 = 9 2 , 
where 9 is the angle between vector k and the z-axis; 9 2 is the mean square scattering 
angle per unit path. The mean square scattering angle ( 9 2 z) at depth z will further be 
denoted by 9 2 , and its square root — by 9 Z . 

We shall also point out that for large values of 9, the solution of (13 . 1 6[) will generally 
be proportional to 1/9 4 [201 [21]. Thus, approximate equation ( 13 .20 j) reflects the fact that 
the large-angle deflections are neglected in a single scattering event. Indeed, for 9 ^> 9 Z , 
the solution (I3.2ip decreases exponentially, while the subsequent term that corresponds 
to the solution of the original equation (13.161) diminishes according to the power law 
9~ A . That is why the angular distribution (13.231) does not describe the characteristics of 
particles scattered at large angles, 9 3> 9 Z [20], [21] . 

The expression for the density matrix in the zeroth approximation ( I3.2ip enables 
one to obtain the correction p^\k,z) in the first-order perturbation theory: 



p 



N 



o 



G{k - k"; 



z - 



<cP i0 \k",z> 



±i 

NC 



a 



1 



NC, 



{(Sn") 2 ,p^Ck",z 



')} 



nucl (0) \{Sn") 2 ,pV ) {k",z') 



G(k - k"; 



z — z 



dQ^,dz' 
fcoui(k , k )p°\k , z )/n Ucl 

,coul 



+N 

Jo J 

+ /nucl,coul(^, k')p m {k\ z')f: Qul (k', k")) dn^Az' 
rz r 

+N 



where G(k 
1 



7l9 2 



z — z 



cxp 



(k\ k h 



G{k - k"; z - z') 

x /nuci,coui(fc", k')p m {k\ z')f: uchcoul (k', k")dn p dn pf dz', 

G(k - k"; z - z') = 
In the general case, the amplitudes for Coulomb and 
3y formula (12. 6p . 



(3.24) 



k"; z — z') is the Green function of (I3.20p : 

(n — n") 2 

'W\z-z'\ 



nuclear scattering are given 

In the case of high energies, the characteristic angles for elastic scattering of 
deuterons by nuclei are 9 n ~ l/kR d 1 (R d is the deuteron radius). As a result, 
in analysing the polarization of high-energy deuteron beams in transmitted geometry 
(figure [[)) using an axially symmetrical detector, the amplitude of scattering due to 
nuclear interaction, which appears in the expression for the density matrix ( I3.24p . can 
be written in the form: 



/nucl,coul(fc, k') = d{9) + d 1 (9)(Sn ) 



(3.25) 



where no is the unit vector directed along the z-axis, which coincides with the direction 
of motion of the initial beam; d(9) and ^1(6*) are the spin-independent and the spin- 
dependent parts of the nuclear amplitude (which allows for the distortion caused to the 
waves incident onto the nuclei by the Coulomb interaction), respectively; 9 is the angle 
between vectors k and k! . 
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The spin structure of the Coulomb scattering amplitude / coul also has the form 
(I3.25p . As stated above, however, for this amplitude, the term proportional to the 
deuteron spin is small and will be dropped further. 



3.2.2. Tensor Polarization of Deuterons Passing Through Matter. By substituting the 
explicit form of the density matrix (12. 5p and the scattering amplitude (I3.25P into the 
solution (I3.24p . one can obtain the characteristics of the deuteron beam, which are of 
interest to us. 

As a result, the number of particles that are registered by the detector with angular 
width Af2 can be written in the form: 



X(z,t?det) 



1 - (-f ' 



h + 6(^tfdet)io 



2 1 

+ g(Qo^oH 



(3.26) 



According to (I3.26p . the intensity of the beam that reaches the detector depends on its 
possible initial tensor polarization Q . 

The parameters £i, £2; and £3 in the small-angle approximation are defined as 
follows: 



£iO,Aiet) = -Na° NC 



1 — exp 



91 



+2itNz J P(x;tfde t , 91) { 2Re[a( X )d*( X )} + \d{x)? )X&X, 


00 

+AttNz [ P(x;$ de t,~9l) ( Re[a(xK(x)] + Md(xK(x)\ ) xd X , 



&Mdet) = 2nNz / P(x;^de t , 9*)\d 1 { X )\ 2 xd X , 



In 



91 . 



9d9, 



(3.27) 



(3.28) 



where the function P(x', $det, 9 2 z ) is defined as 

p(x',***, 91) = = J - P (-^L 

where Io(y) is the modified Bessel function of the zeroth order. 

Let us discuss now what polarization this beam has. Write the expressions for 
the number of particles that are registered by the detector and have the initial spin 
projections M = and M = ±1, respectively: 



X (0) (Mdet) 
X (±1) (z,tf d et) 



1 

Vdct 



-'0 + ?1 J > 



(0) 



1 — exp 



91 



it i) + (6 + 6 + wt } 



(3.29) 
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For the sake of simplicity, let us consider the case when an unpolarized beam is 
incident onto the target. Using the definition (I3.3p . find tensor polarization of the beam 
registered by the detector: 

2 



Pz 



exp { - 



(6 + 6). 



(3.30) 



Here we have retained the terms proportional to the first power of the small 
quantities £ 2 and 6- 

Using the explicit form of the parameters 6; 6 m (13.27ft . the expression for tensor 
polarization of the beam can be written in the form: 



Pz 



dct J 



cr°)z 



-N 



47T 



1 - exp 
2tt 



P( X ,tfdet, ^ 2 )Re[a(x)^(x)]xdx 



X 



(2Re[d(xK(x)] + K(x)| 2 ) X^X 



(3.31; 



Equation (13.3 ip simplifies appreciably for the two limiting values of the detector 
angle: i?det "C # 2 and ^det ^> 9 Z . As has been stated above, the solution (I3.2ip does 
not describe the properties of the beam that is scattered at angles 9 > 9 Z . But the 
integral characteristics of deuterons that are obtained from this solution can be used for 
i^det > 9 Z . The influence of particles single-scattered due to the Coulomb interaction 
can be neglected because the major contribution to (I3.28P comes from the small values 
of 9. 

When $det is much less than 9 Z , the expression for p zz can be written in the form: 
2 
3' 



Pzz(z, ^det < 9 Z ) 



-N(a. 



±i 

tot 



°tot) z 



+^nNz J exp (-== j Re[a(x)dl(x)]x^X 



4 

+ -ttNz 
3 



J exp (2Re[d(x)dt( X )} + Mi(x)| 2 ) xd X - 



(3.32) 



According to (13. 32 p . as a result of multiple scattering, tensor polarization becomes 
independent of $dct for small angles of the detector. Equation ( I3.32p contains a 
scattering-related contribution (the second and third terms), which appears in (I3.32p as 
a result of deuteron rescattering from the direction of k ^ k into the direction of vector 

/co- 
in the other limiting case, when $d et is much greater than 9 Z , the terms related to 
scattering in (13 .311) are equal to the scattering cross sections, and so the expression 
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between the square brackets vanishes. Tensor polarization is determined by the 
difference between the cross sections af 1 and of alone, as it occurs in the case of a 
thin target. 

Let us represent the expression (I3.3ip for tensor polarization in the form: 

M^tfdet) = -\N{a? - cr°)z - ^N(af - a el )zR(z,# det ), (3.33) 

where the function H is defined as 

1 4tt 



H(z,tf dct ) = l- 



1 - exp ( -^f 



P(x,#det, 9^)Re[a(x)dl(x)}xdx 



\0 



P(x,^ d et, 61) (Re[d(xK( X )} + \di(x)\ 2 /2) X d X ) • (3.34) 



o 



Note that when passing to the limit of a thin target, i.e., z — > 0, one should 
take into account that in this case, the function P(x;$det, $ 2 ) turns to P(x; $det, 0) = 



°dct 

J 5(i3 — x)di9, where 6 is the delta function. 



According to (13. 33 p . the dependence of tensor polarization of a beam on the target 
thickness and the angle $det is determined by the function H(z, ^det)- 

Let us recast the equality for the mean square angle of deuteron scattering at depth 
z in the form 1191: 



9 2 z = 16nZ 2 {^\ \n{137 Z-V^^j-AC, (3.35) 

where AC is the target thickness in g/cm 2 , Na is the Avogadro's number. 

The explicit form of the d dependence on 9 in the range of small scattering angles 
for the case of structural particles was derived, for example, in [22]. To estimate the 
main parameters, let us consider the simplest form of the d and d\ dependence on 9, 
namely, 

d{9) = d(0) exp (-A; 2 P^ 2 /4) , d x {9) = d x (0) exp (-A; 2 P^ 2 /4) ,(3.36) 

where Pd is the deuteron radius. 

We also used the fact that according to the eikonal- approximation calculations [23] , 
in the range of energies of incident deuterons from 400 to 800 MeV, for the amplitude 
d we have Imd = 0.75 • 10~ n cm, Red = —0.6 • 10 _12 cm, and the imaginary part of the 
amplitude d\ relates to its real part as Imrfi/Re^i ~ —10, with the spinless part of the 
nuclear amplitude / nuc i being much larger than its spin part. 

For the considered energies of deuterons scattered by a screened Coulomb potential, 
the first Born approximation to the amplitude is satisfactory to use as a real part of 
the Coulomb amplitude a(9), and the second term of the Born series can be used as its 
imaginary part accordingly. 

Figure[2]and figure [3] present the results of calculations of the function H for different 
values of f^det and different deuteron energies when a beam is scattered by a carbon filter. 
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Figure 2. Function H versus the detector angle for deuteron energy of 500 MeV. Solid 
curve corresponds to the carbon target thickness ACi = 2 g/cm 2 , dashed curve - to 
AC 2 = 10 g/cm 2 , dotted curve — to AC 3 = 20 g/cm 2 . 
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Figure 3. Function H versus the deuteron energy for the detector angle t?dct = 5 • 10~ 3 
rad. Solid curve corresponds to the carbon target thickness ACi = 2 g/cm 2 , dashed 
curve — to AC2 = 10 g/cm 2 , dotted curve — to AC3 = 20 g/cm 2 . 

The results of calculations given in figure [2] demonstrate that owing to multiple 
Coulomb scattering, the value of H is less than unity, H < 1, even in a narrow-angle 
geometry (tfdet <C 6 Z ). Let us recall that this occurs because of the rescattering in 
the direction of vector k of particles that have been scattered in different directions. 
For AC = 2, 10, and 20 g/cm 2 , the values of 6 Z = \fWz are 6 Z (AC 1 ) = 4.5 mrad, 
9 Z (AC 2 ) = 10 mrad, and 9 Z (AC 3 ) = 14 mrad, respectively. As a result, in this 
limiting case, the values of H for deuterons registered by the detector with angular 
width $ det < 6 Z for target thicknesses AC = 2, 10, and 20 g/cm 2 are H(AC X ) = 0.68, 
H(AC2) = 0.59, and H(AC3) = 0.54, respectively. As is seen, these values differ from 
unity by almost a factor of two, and this difference is the largest for the thickest target. 

With growing angle &det, the particles scattered at the angles 6 < t^et gain in 




T 




AC 1 =2g/cm 2 . 

AC 2 =10g/cm 2 • 

AC 3 =20g/cm 2 ; 
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significance, which leads to a decrease of H. In a wide-angle geometry, which for the 
fast-particle scattering is defined from the condition $ det 3> 9 e g (0 e g = 6 Z when 6 Z > 9 n , 
and # e ff = 6n when 9 Z < 6 n ), we have H — > 0. As a result, tensor polarization in this 
limiting case is determined by the difference between the inelastic parts of the total 
scattering cross sections for the states with M = ±1 and M — 0. 

Now, let us consider the behaviour of the function H for different values of energy at 
a fixed value of fi^et (choose $det to be $det = 5 ■ 1CT 3 rad). According to figure [31 when 
the deuteron energy grows, the value of H first increases and then starts to diminish. 

It should be emphasized that with due account of multiple Coulomb scattering, 
tensor polarization is a nonlinear function of the target thickness z ( I3.3ip . For large 
values of $det; however, p zz becomes a linear function of z (13. 10p . With increasing 6 Z , 
expression (I3.32H approaches the expression for tensor polarization ( 13.101) . This means 
that for large values of z or for low energies of the incident deuterons, the dependence 
of the tensor polarization on the detector's angle becomes insignificant (it is necessary 
that the parameter 9 Z should be much greater than the characteristic angle 9 n of nuclear 
scattering). In this case, tensor polarization can be considered linear in z for any values 

of $det- 

4. Conclusion 

It has been shown that the magnitude of tensor polarization of the deuteron beam, which 
arises from the spin dichroism effect, depends appreciably on the angular width i?det of 
the detector that registers the deuterons transmitted through the target. Even when 
the angular width of the detector is much less than the mean square angle of multiple 
Coulomb scattering, the beam's tensor polarization depends noticeably on rescattering. 
In the case when the angle d^et is much larger than the mean square angle of multiple 
Coulomb scattering (as well as than the characteristic angle of elastic nuclear scattering), 
tensor polarization is determined by the difference between the total reaction cross 
sections for deuteron-nucleus interaction, af 1 — of- Elastic scattering processes here 
make no contribution to tensor polarization. 
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